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Abstract. The notion of p-adic multiresolution analysis (MRA) is intro- 
duced. We discuss a "natural" refinement equation whose solution (a refin- 
able function) is the characteristic function of the unit disc. This equation 
reflects the fact that the characteristic function of the unit disc is a sum 
of p characteristic functions of mutually disjoint discs of radius p~^ ■ This 
refinement equation generates a MRA. The case p = 2 is studied in de- 
tail. Our MRA is a 2-adic analog of the real Haar MRA. But in contrast 
to the real setting, the refinable function generating our Haar MRA is 1- 
periodic, which never holds for real refinable functions. This fact implies 
that there exist infinity many different 2-adic orthonormal wavelet bases in 
'C^(Q2) generated by the same Haar MRA. All of these bases are described. 
We also constructed multidimensional 2-adic Haar orthonormal bases for 
'C^(Q2) by means of the tensor product of one-dimensional MRAs. A cri- 
terion for a multidimensional p-adic wavelet to be an eigenfunction for a 
pseudo-differential operator is derived. We proved also that these wavelets 
are eigenfunctions of the Taibleson multidimensional fractional operator. 
These facts create the necessary prerequisites for intensive using our bases 
in applications. 



1. Introduction 

According to the well-known Ostrovsky theorem, any nontrivial valuation 
on the field Q is equivalent either to the real valuation \ ■ \ or to one of the p- 
adic valuations \ ■ \p. This p-adic norm | ■ |p is defined as follows: if an arbitrary 
rational number a; 7^ is represented as x = p'^^, where 7 = 7(x) G Z and 
the integers m, n are not divisible by p, then 

(1.1) \x\p=p-\ x^O, \0\p = 0. 

The norm | ■ \p satisfies the strong triangle inequality + ^ max(|x|p, \y\p) 
and is non-Archimedean. The field Qp of p-adic numbers is defined as the 
completion of the field of rational numbers Q with respect to the norm | ■ L. 
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Thus there are two equal in rights universes: the real universe and the p- 
adic one. The latter has a specific and unusual properties. Nevertheless, there 
are a lot of papers where different applications of p-adic analysis to physi- 
cal problems, stochastics, cognitive sciences and psychology are studied [6]- 
[To] . [13] [19] . [35] [37] (see also the references therein). In view of the 
Ostrovsky theorem, such investigations are not only of great interest in itself, 
but lead to applications and better understanding of similar problems in usual 
mathematical physics. 

Since there exists a p-adic analysis connected with the mapping Qp into 
Qp and an analysis connected with the mapping Qp into the field of complex 
numbers C, one considers two corresponding types of p-adic physical models. 
For the p-adic analysis related to the mapping Qp — C, the operation of 
differentiation is not defined, and as a result, large number of models connected 
with p-adic differential equations use pseudo-differential operators (see the 
above-mentioned papers and books). In particular, fractional operator are 
extensively used in applications. A very important fact that the eigenf unctions 
of a one- dimensional fractional operator D"' form an orthonormal basis for 
£^(Qp) was observed by V.S. Vladimirov, I.V. Volovich, E.I. Zelenov (see [35]). 
S. V. Kozyrev [20j found an orthonormal compactly supported p-adic wavelet 
basis for £^(Qp): 

(1.2) ek.,,aix) = p-^'\p{p'^k{pix - a))^l{\pix - a|p), x E Qp, 

k = 1,2, ... ,p — 1, j E Zi, a E Ip = Qp/Zp. Wavelets (11.21) are also eigenfunc- 
tions of the one-dimensional fractional operator D"': 

D'^Ok-jaix) = p"(^-^)4;,.(x), a; G Qp, a E C. 

Some wavelet- type systems generahzing (11.21) were suggested by S. V. Kozyrev 
[H], [22], A. Yu. Khrennikov and S. V. Kozyrev [16j, [17\, J. J. Benedetto 
and R. L. Benedetto [8j, R. L. Benedetto [9]. Multidimensional p-adic bases 
obtained by direct multiplying out the Kozyrev's wavelets (11.21) were considered 
in [3]. The authors of [18] found the following new type of p-adic wavelet basis: 

(1.3) 9^:;;l{x)=p-^/\p{s{p^x-a))n{\p^x-a\p), xEQp, 

where m > 1 is a fixed positive integer; s = p~"^{so + sip + ■ • • + Sm-ip"*"^), 
Sr = 0, 1, . . . ,p — 1, r = 0, 1, . . . , m — 1, So 7^ 0; j E Z, a E Ip. It turned 
out that wavelets (II. 2p and their generahzations are eigenfunctions of |j-adic 
pseudo-differential operators [3]- [5], [16], [17], [18], [20] - [22]- Moreover, a 
necessary and sufficient conditions for a class of p-adic pseudo-differential op- 
erators (12.151) (including fractional operator (12.231) ) to have wavelets (II. 2p and 
(II. 3p as eigenfunctions was derived in ||3j , [ISJ . So, wavelets play an important 
role for application of p-adic analysis and gives a new powerful technique for 
solving p-adic problems. Nevertheless, in the cited papers, there was no any 
attempt to create a theory describing common properties of wavelet bases and 
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giving methods for their finding. The goal of this paper is to start development 
of such a theory. 

It's interesting to compare appearing first wavelets in p-adic analysis with 
the history of the wavelet theory in real analysis. In 1910 Haar pL2j constructed 
an orthogonal basis for £^(M) consisting of the dyadic shifts and scales of one 
piecewise constant function. A lot of mathematicians actively studied Haar 
basis, different kinds of generalizations were introduced, but during almost the 
whole century nobody could find another wavelet function (a function whose 
shifts and scales form an orthogonal basis). Only in early nineties a general 
scheme for construction of wavelet functions was developed. This scheme is 
based on the notion of multiresolution analysis (MRA in the sequel) introduced 
by Y. Meyer and S. Mallat p8J, [26j, [27j|. Smooth compactly supported wavelet 
functions were found in this way, which has been very important for various 
engineering applications. In the present paper we introduce MRA in £^(Qp) 
and study a concrete MRA for p = 2 being an analog of Haar MRA in £^(]R). 
The same scheme as in the real setting leads to a Haar basis. It turned out 
that this Haar basis coincides with Kozyrev's wavelet system (11. 2p . However, 
2-adic Haar MRA is not an identical copy of its real analog. We proved that, 
in contrast to Haar MRA in £^(M), there exist infinity many different Haar 
orthogonal bases for C?{'^2) generated by the same MRA. 

The paper is organized as follows. 

In Sec. [21 we recall some facts from the p-adics. The basic results on the the- 
ory of the Bruhat-Schwartz distributions are given in Subsec. 12.11 (see [11] , [33] , 
[M] , |35j). Some facts from the theory of the p-adic space ^'(Qp) of Lizorkin 
distributions [3] are summarized in Subsec. 12. 2[ In Subsec. I2.3[ \2A\ we re- 
call some facts [3] related to the multidimensional pseudo-differential opera- 
tors (I2.15P in the space of Lizorkin distributions ^'(Qp). In particular, the 
multidimensional fractional operators introduced by Taibleson [33l §2], [3ll 
III. 4.] are discussed. The spaces $'(Qp) is a "natural" domain for the class 
of pseudo-differential operators (I2.15p . The spaces $'(Qp) are invariant un- 
der our pseudo-differential operators. It is appropriate to mention here that 
the class of our operators includes the pseudo-differential operators studied 

in iM], m- 

In Sec. [31 a notion of p-adic MRA is introduced (Definition 13.11) . In Sub- 
sec. 13.21 we discuss refinement equation (13.70 : 

1 r 

0W = y]0( -a; - -), a; e Qp, 

whose solution (a refinahle function) (p is the characteristic function i7(|a;|p) 
of the unit disc. The conjecture to use the above equation as a refinement 
equation was proposed in [12]. The above refinement equation is natural and 
reflects the fact that the characteristic function of the unit disc Bq = {x : 
\x\p < 1} is represented as a sum of p characteristic functions of the mutually 
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disjoint discs 5„i(r) = {x : || - < l}, r = 0, 1, . . . ,p - 1 (see (lO) ). 
In Subsec. 13.31 the 2-adic Haar MRA is constructed. Namely, we proved that 
the refinable function (j){x) = il(^\x\2) generates a MRA, which is an analog of 
the classical Haar MRA. It is shown that a 2-adic analog of the real wavelet 
function generated by Haar MRA generates an orthonormal basis fl3.14p for 
^^(Q2)- This basis coincides with Kozyrev's one (11. 2p for p = 2. We proved 
that Kozyrev's basis is not a unique orthonormal wavelet basis generated by 
2-adic Haar MRA. 

In Sec. m infinity many 2-adic Haar wavelet functions generating different 
orthonormal bases for £^(Q2) are found. Explicit formulas for Haar wavelet 
functions s eN, (Theorem 14. 11) are given in Subsec. 14. 1[ All real wavelet 
functions are described for s = 1, 2 in Subsec. 14.21 

In Sec. El we study multivariate Haar bases. A general scheme for construc- 
tion of 2-adic multidimensional separable MRA is described in Subsec. 15.11 . 
According to this scheme, separable 2-adic Haar wavelets (15.71) in ^^(Qg) are 
constructed in Subsec. 15. 2[ 

Sec. |6] is devoted to the spectral theory of pseudo-differential operators 
(12.151) . Theorems 16.11 16.21 state a criterion for pseudo-differential operators 
(12.151) to have multidimensional wavelets (15. 5p or (15.71) as eigenfunctions. In 
particular, according to Corollaries 16.11 16.21 the multidimensional wavelets 
(15. 5p and (15. 7p are eigenfunctions of the multidimensional fractional operator 



2. Preliminary results in ^-adics 

2.1. p-Adic distributions. Here and in what follows, we shall systematically 
use the notations and the results from [35] and [lit Ch.II]. Let N, Z, C be the 
sets of positive integers, integers, complex numbers, respectively. No := {0}UN. 

The field of p-adic numbers is denoted by Qp. The canonical form of any 
p-adic number a; 7^ is 

(2.1) X = p'^{xo + XiP + X2p'^ -\ ), 

where 7 = 7(x) G Z, Xj = 0, 1, . . . ,p — 1, xq 7^ 0, j = 0, 1, . . . . The series 
is convergent in the p-adic norm (II. ip . and one has \x\p = p~^. By means of 
representation (12. ip . the fractional part {x}p of a number a; G Qp is defined as 
follows 

l^-^; l^ip \ pT(xo + a;ip + a;2P=^ + --- + a;|^l_ipl^l-i), if 7(x) < 0. 

The function 
(2.3) xpi^x) = e^-'^^^^- 

for every fixed ^ G Qp is an additive character of the field Qp, where {-jp is a 
fractional part (12.20 . 
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The space Qp := Qp x • ■ ■ x Qp consists of points x = (xi, . . . , Xn), where 
Xj e Qp, j = 1,2 . . . ,n, n > 2. The p-adic norm on Qp is 

(2.4) \x\p = max |x,-L, x e 

l<j<n ^ 

where I defined by ffTTD . 

Denote by Bl!j{a) = {x G Qp : |x — a|p < p^} the ball of radius with the 
center at a point a = (ai, . . . , a„) G Qp and by 5'"(a) = {x G Qp : |x — a|p = 
pT| = B!^{a) \ i?"_]^(a) its boundary (sphere), 7 G Z. For a = 0, we set 
-B"(0) = and 5'"(0) = 5". For the case n = 1, we will omit the upper index 
n. It is clear that 

(2.5) B:;{a)=B^{a^)x---xB^{an), 

where B^{aj) = {xj : \xj — aj\p < p'^} C Qp is the disc (one-dimensional ball) 
of radius p^ with the center at a point G Qp, j = 1, 2 . . . , n. 

Any two balls in Qp either are disjoint or one contains the other. Every 
point of a ball is its center. 

According to [35l 1. 3, Examples 1,2.], a one-dimensional disc B^ is repre- 
sented as a sum of p'^'^ mutually disjoint discs By (a), 7' < 7: 

(2.6) B^ = ByU (^UaBy{a)^, 

where a = and a = a_rP~''' + cb~r+iP~^'^^ + ■ — I- a^y^ip^''''^^ are the centers 
of the discs By (a), r = 7, 7 — 1, 7 — 2, . . . , 7' + 1, < aj < p — 1, a^r 7^ 0. 
In particular, the disc Bq is represented as a sum of p mutually disjoint discs 

(2.7) i?o = i?-iU (u^li^i?_i(r)), 

where B_i{r) = {x E So : Xq = r} = r + pZp, r = 1, ... ,p — 1; B_i = 
{\x\p < P~^} = P'^p] and So = {|x|p = 1} = U^Z]^-B_i(r). Here all the discs are 
mutually disjoint. Coverings (12. 6p and (12. 7p are called the canonical covering 
of the discs Bo and B^, respectively. 

There exists the Haar measure dx on Qp. This measure is positive, invariant 
under the shifts, i.e., d{x + a) = dx, and normalized by Jj^| dx = 1. The 
invariant measure dx on the field Qp is extended to an invariant measure 
d^x = dxi ■ ■ ■ dxn on Qp in the standard way. 

A complex- valued function / defined on Qp is called locally- constant if for 
any x G Qp there exists an integer /(x) G Z such that 

f{x + y) = f{x), yeB^^y 

Let i^(Qp) and '^'(Qp) be the linear spaces of locally-constant C-valued func- 
tions on Qp and locally-constant C-valued functions with compact supports 
(so-called test functions), respectively [351 VI. 1., 2.]. If y9 G ^'(Qp), according 
to Lemma 1 from [SS', VI. 1.], there exists / G Z, such that 

^{x + y) = v{x), yeB^, x G Q" 
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The largest of such numbers I = l{ip) is called the parameter of constancy of 
the function (p. Let us denote by Pjy(Qp) the finite-dimensional space of test 
functions from I?(Qp) having supports in the ball and with parameters 
of constancy > / [35i, VI. 2.]. The following embedding holds: PjylQp) <^ 
Pj;,(Qp, N <N',l> I'. Thus V{Q;) = limind^^oolimindz^_ooI^5v(Qp)- 
The space V{Qp) is a complete locally convex vector space. According to [Ssl 
VI, (5. 2')], any function ip G r'jy(Qp) is represented in the following form 

pn{N-l) 

(2.8) ipix)= J2 ^(Ofi(P"'|a;-c'^|,), xgQ^, 

where fl{p~''\x — c'^\p) are the characteristic functions of the mutually disjoint 
balls Bi{c'^), and the points = (c^', ■ ■ ■c'Q G do not depend on (p. 

Denote by T>'{Qp) the set of all linear functional (p-adic distributions) on 

viQ;) m VI.3.]. 

Let us introduce in V{Qp) a canonical S-sequence 6k{x) = p"'^Q{p''\x\p), and 
a canonical 1-sequence Q{p~''\x\p) (the characteristic function of the ball B^), 
k e Z, X e Qp, where 

'1, < t < 1, 
0, t > 1, t G M. 



(2.9) n{t) 



It is clear that 6k ^ S, k ^ oo in V'{Qp), where 6 is is the Dirac delta function, 
and Q{p-^\ ■ \p) ^ 1, k ^ oo in ^(Qp (see [351 VI.3., VII.l.]). 
The Fourier transform of G V{Qp) is defined by the formula 



F[m)= [ Xp{^-xMx)rx, 



where Xp{^ ' ^) = Xpi^i^i) ' ' 'XpHnXn) = e^'^'^'^^=^^^^^^^'^] ^ ■ a; is the scalar 
product of vectors and Xpi^j^j) are additive characters (12. 3p . The Fourier 
transform is a linear isomorphism T>{Qp) into V{Qp). Moreover, according 
to [331 Lemma A.], [Ml 111,(3.2)], [351 Vn.2.], 

(2.10) ^{x)evUQ;) iff F[^{x)]{OeVz!'{Q;). 

The Fourier transform F[f] of a distribution / G 2^'(Qp) is defined by 

(F[/],^) = (/,FM), V^GP(Q^). 

Let y4 be a matrix and b G Q^. Then for a distribution / G T>'{Qp) the 
following relation holds [351 VII, (3.3)]: 

(2.11) F[f{Ax + 6)](0 = I det A\;'xp{ - A'^b ■ ^)F[f{x)] {A-'^), 
where det A 7^ 0. According to p5l, IV, (3.1)], 

(2.12) F[Q{p~''\-\p)]{x) = 6k{x), keZ, xeQ;. 
In particular, F[Q{\ ■ |p)](x) = i7(|a;|p). 
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2.2. The p-adic Lizorkin spaces. Let us introduce a space of the p-adic 
Lizorkin test functions (see |3], [1]) 

$(Q;) = {0:0 = i^M, ^G^(Q;)}, 

where ^(Q^) = {^(0 G ^^(Q^) : ^(0) = 0}. Here ^(Qp,<l>(Q^) C V{Qp. 
The space $(Qp) can be equipped with the topology of the space V{Qp) which 
makes $ a complete space. 
In view of f l2^ . 

(2.13) e m;) ^ [ d'^x = 0, e v{q;). 

Moreover, G ^'^(Q;^) H $(Qp, i.e., 0(x) = 0, iff ^ = ^ 

vz^iQ;) n v[/(Qp, i.e., ^(0 = 0, e e 5^^". 

Let $'(Qp) and ^'(Qp) denote the topological dual of the spaces $(Qp) and 
\E'(Qp), respectively. We call $'(Qp) the space of p-adic Lizorkin distributions. 

If and are subspaces of functionals in T>'{Q^) orthogonal to ^&(Qp) 
and $(Q;J), respectively, then = {/ e P'(Qp : f = C5, C E C} and 
$^ = {/ e P'(Qp") : / = C, C G C}. According to [3|, $'(Qp = V'{Q;)/<^>^ 
and ^'(Qp) = P'(Qp/^-^. Thus the space $'(Qp) can be obtained from 
by "sifting out" constants. 

We define the Fourier transform of / G $'(Qp) and g G ^E''(Qp) respectively 

by 

(F[/],^) = (/,F[V^]), VV^GV^(Qp, 

Since F[$(Qp] = ^(Q;j) and F[^(Q;j)] = $(Q;j), formulas fimj) give well de- 
fined objects. It is clear that F[$'(Qp] = ^'(Qp and F[^'{Q^)] = $'(Qp [3]. 

Recall that the usual Lizorkin spaces were studied in the excellent papers of 
P. I. Lizorkin [2l], [25| (see also ^Slj, [32]). 

2.3. Pseudo-differential operators in the Lizorkin space. Consider the 
following class of pseudo-differential operators A in the Lizorkin space of the 
test functions $(Qp) defined by 

{A<p){x) = F-'[AF[<p]]{x) 

(2.15) = / / Xp{{y-x)-^)Am{y)d''^d^y, 0e$(Qp, 

with symbols A G S{Qp \ {0}), which were introduced in [3]. 

Lemma 2.1. ( [3J) The Lizorkin space $(Qp) is invariant under the pseudo- 
differential operators fl2.15p . Moreover, A($(Qp)) = $(Qp). 
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Given pseudo-differential A with a symbols A, define the conjugate operator 
A'^ by 

(2.16) {A^<P){x)=F-'[Ai-)Fmx)= [ Xp{-x ■ OA-OmiO d""^ 

and introduce the corresponding operator A in the space of Lizorkin distribu- 
tions: if / e $'(Qp), then 

(2.17) (A/,0) = (/,A^</>), y<f>em;)- 

It is clear that 

(2.18) Af = F~'[AF[f]]e^'iQ;), 

i.e., the space of Lizorkin distributions $'(Qp) is invariant under pseudo- 
differential operators A. Moreover, by Lemma [2?n A($'(Qp)) = $'(Qp). 

If A,B are pseudo-differential operators with symbols A,B& ^{Qp \ {0}) 
respectively, then the operator AB is well defined and represented by the 
formula 

{AB)f = F-'[ABF[f]]e<l>\Q;). 

If v4(^) 7^ for all ^ G Qp \ {0}, then operator we the pseudo-differential 
operator 

A-'f = F-'[A-'F[f]], /G<I>'(Q;). 

is, evidently, the inverse operator to A. 

Thus the family of pseudo-differential operators A whose symbols A do not 
vanish on \ {0}, forms an Abelian group. 

2.4. Multidimensional fractional operator. Let us consider a pseudo- dif- 
ferential operator D°' with the symbol A{^) = \C,\p, a G C, a 7^ —n. By 
definition fl2J[5|) . 

(2.19) (D»(x) = F-i[|-|p^M(-)](^), 0e$(Qp"). 

These operators were introduced by Taibleson [33l §2], [Ml III. 4.] as operators 
in the space of distributions I^'(Qp)- The relation (12.191) can be rewritten 

as the following convolution (see [35l VII. 1.]): (D'^0)(x) K_a(x) * 0(x) = 
(fi:_„(0, 0(a; - 0). xeQ^,^e $(Qp), where the distribution (from V'Iq^)) 

(2.20) K^{x) = ^T^, a^O, n, xe 

is called the multidimensional Riesz kernel, r^\a) is the n-dimensional F- 
function (see for details, see [33], [Ml III], [35l VIII]). 

The Riesz kernel has a removable singularity at a = and according to [33l 
§2], [Ml III.4.], [351 VIII.2], we have 

(2.21) Ko{x) lim Ka{x) = S{x). 
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According to |3], [1], one can similarly introduce a Lizorkin distribution Kn{-) 
by 

(2.22) Kn{x) =^ lim Ka{x) = -i— ^ — log \x\p. 

a^n lOgp 

Due to fl2.2ip , fl2.22p , we can define the operator fl2.19p in the Lizorkin space 
of test functions (f) G $(Qp) by 

(2.23) (D»(x)t//€_,(x)*0(x) = (/€_,(O,0(x-O), xeQ;, aeC. 

According to Lemma 12.11 the Lizorkin space $(Qp) is invariant under the 
Taibleson fractional operator and L)"($(Qp) = $(Q;j). By (HJM . (1^371) . 
(L)")T = and for any / G $'(Qp) we have 

(2.24) t/ V0 G $ (Q;). 

It is clear that D"(<l''(Q;j)) = $'(Qp. 

It is easy to prove that Ka(x) * fi;^(x) = Ka+f^^x), a, /? G C, holds in the sense 
of the Lizorkin space $'(Qp) (see [3J). Consequently, the family of operators 
D", a G C, in the Lizorkin space forms an Abelian group: if / G $'(Qp) then 
/ = D^D"/ = = a,/5GC. 

3. MULTIRESOLUTION ANALYSIS (ONE-DIMENSIONAL CASE) 

3.1. p-Adic multiresolution analysis. Consider the set 

Ip = {a = p^'' (flo + CLiP H h a^-ip^'^) : 

(3.1) 7 e N; a,- = 0, 1, . . . ,p - 1; J = 0, 1, . . . , 7 - 1}- 

This set can be identified with the factor group Qp/Zp. 

It is well known that Qp = BqU W^^-^S^, where 5*^ = {x G Qp : \x\p = p'^}. 
Due to (12. ip . X G S"^, 7 > 1, if and only if x = X-^p~^ + X-^+ip^'''^^ + ■ ■ ■ + 
X-ip~^ + ^, where 7^ 0, ^ G Bq. Since X-^p"'^ + X-^+ip~'^~^^ + ■ ■ ■ + X-ip~^ G 
Ip, we have a "natural" decomposition of Qp to a union of mutually disjoint 
discs: 

Qp = U 5o(a). 

So, Ip is a "natural" group of shifts for Qp, which will be used in the sequel. 

Definition 3.1. A collection of closed spaces Vj C £^(Qp), j G Z, is called a 
multiresolution analysis (MRA) in £^(Qp) if the following axioms hold 

(a) Vj C Vj+i for all j G Z; 

(b) U Vj is dense in C^iQp); 

(c) Q/^ = {oy, 

(d) /(■) G V,- ^ /(p-i-) G for all j G Z; 
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(e) there exists a function (p E Vq such that the system {0(- — a), a G Ip} is 
an orthonormal basis for Vq. 

The function (p from axiom (e) is called refinable or scaling. It follows im- 
mediately from axioms (d) and (e) that the functions p^^'^cpip"^ ■ —a), a E Ip, 
form an orthonormal basis for Vj, j E Z. 

According to the standard scheme (see, e.g., [301 §1-3]) for construction of 
MRA-based wavelets, for each j, we define a space Wj {wavelet space) as the 
orthogonal complement of Vj in V^+i, i.e., 

(3.2) = V, © W,, J E Z, 
where Wj -L Vj, j E Z. It is not difficult to see that 

(3.3) f eW,^ f{p-'-)EWj+i, for all jEZ 

and Wj _L Wk, j 7^ k. Taking into account axioms (b) and (c), we obtain 

(3.4) Wj = C^iQp) (orthogonal direct sum). 

If now we find a function ip E Wq such that the system {ip^x — a),a E Ip} 
is an orthonormal basis for Wq, then, due to (13.31) and (13. 4p . the system 
|p772^(p-i . —a), a E Ip,j E Z}, is an orthonormal basis for £^(Qp). Such 
a function is called a wavelet function and the basis is a wavelet basis. 

3.2. p-Adic refinement equation. Let be a refinable function for a MRA. 

As was mentioned above, the system {p^/'^(p(p^^ ■ ~a),a E Ip}, is a basis for 
Vi. It follows from axoim (a) that 

(3.5) (p =^aa(pip~^ ■ -a), aa E C. 

We see that the function is a solution of a special kind of functional equation. 
Such equations are called refinement equations. Investigation of refinement 
equations and their solutions is the most difficult part of the wavelet theory in 
real analysis. 

A natural way for construction of a MRA (see, e.g., [501 §1.2]) is the fol- 
lowing. We start with an appropriate function (p whose integer shifts form 
an orthonormal system and set Vj = span{0(p~-5' ■ —a) : a E Ip}, jEZ. It 
is clear that axioms (d) and (e) of Definition 13.11 are fulfilled. Of course, 
not any such a function (p provides axiom (a). In the real setting, the re- 
lation Vq C Vi holds if and only if the refinable function satisfies a refine- 
ment equation. Situation is different in p-adics. Generally speaking, a re- 
finement equation (13. 5p does not imply the including property Vq C Vi. In- 
deed, we need all the functions (p{- — h), b E Ip, to belong to the space Vi, 
i.e., the identities (p{x — b) = J2aeip (^a,b(p{p~^^ ~ should be fulfilled for 
all b E Ip. Since p~^b + a is not in Ip in general, we can not state that 
(p{x — b) = Xlae/ '^a,b0(P~^3; — p^^b — a) E Vi for all b E Ip. Nevertheless, 
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some refinable equations imply including imply property, which may happen 
because of different causes. 

The refinement equation refiects some "self- similarity". The structure of the 
space Qp has a natural "self-similarity" property which is given by formulas 
(12. 6p . (12. 7p . By (12.71) . the characteristic function of the unit disc Bq 

is represented as a sum of p characteristic functions of the mutually disjoint 
discs _B_i(r), r = 0, 1, . . . ,p — 1, i.e., 



(3.6) r2(|s|p) = r2(p|a; — r|p) = -a; a 

r=0 r=0 ^ ^ ^ 

Thus, in p-adics, we have a natural refinement equation (13.51) : 



P-^ 1 

(3.7) (f)(x) = (pf-x xe 



whose solution is 0(x) = This equation is an analog of the refinement 

equation generating the Haar MRA in real analysis. 

3.3. Construction of 2-adic Haar multiresolution analysis. Using the 
refinement equation (13. 7p for p = 2 

(3.8) 0(a;) = 0(^ia;)+0(^la;_l), a; e Q2, 

and its solution, the refinable function 0(x) = fi(|x|2), we construct a 2-adic 
multiresolution analysis. 
Set 



(3.9) Vj = span{(j){2-^x-a) -.aeh}, J eZ. 

It is clear that axioms (d) and (e) of Definition 13.11 are fulfilled and the system 
{2-^/^0(2~'' ■ —a), a G Ip} is an orthonormal basis for V^-, j E Z. Since the 
numbers 2~^b, 2~^b + 2~^ are in I2 for all b G I2, it follows from the refinement 
equation (13.81) that Vq C Vi. By the definition (13. 9p of the spaces Vj, this yields 
axiom (a). Due to the refinement equation (13.80 . we obtain that Vj C V^+i, 
i.e., the axiom (a) from Definition 13. II holds. 

Note that the characteristic function of the unit disc r2(|x|2) has a wonderful 
feature: f2(| ■ +^|2) = ^(| ■ I2), for all ^ G Z2 because the p-adic norm is non- 
Archimedean. In particular, Q{\ ■ ±1|2) = fi(| ■ I2), i.e., 

(3.10) 0(x± 1) = 0(a;), Vx G Q2- 
Thus is a 1-periodic function. 



Proposition 3.1. The axiom (b) of Definition [3TT] holds, i.e., Uje^V^- 
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Proof. According to (12.81) . any function ip G 1^(02) belongs to one of the spaces 
^jvlQs); and consequently, is represented in the form 

2N-1 

(3.11) ^{x) = J2^ic'')n{2-'\x-c''\2), xgQs, 

1>=1 

where e Bn, z/ = 1,2,... 2^"'; / = l{ip), N = N{ip); I G Z. Taking 
into account that n{2-^\x - c'^y = n{\2^x - 2^d'\2) = i{2^x - 2^d'), we can 
rewrite (13.111) as 

2JV-i 

(^(x) = ^ a^0(2'a; - 2^0") x G Qa, c" G 5^, G C. 

Since any number 2'c'' can be represented in the form 2^c'^ = a'^ + b'^, G /a, 
b" E Z2, using (I3.10p . we have 

2N-1 

(p{x) = ^ a,y(f){2''x — a"), x G Q2, G /a, G C, 
i/=i 

i.e., (p{x) G Thus any test function (p belongs to one of the space Vj, 

where j = ji(p), j G Z. 

Since the space ^^(Qa) is dense in ^^(Qa) [35, VI.2], approximating any 
function from ^^(Qa) by test functions ip G V{Q2), we prove our assertion. □ 

Proposition 3.2. The axiom (c) of Definition \3.1\ holds, i.e., (IjfzzVj = {0}. 

Proof. Assume that ^jfLiYj 7^ {0}. Then there exists a function / G ^^(Qa) 
such that 11/11 7^ and / G V^- for all j G Z. Hence, due to (13.91) . we have 
f{'^) = Eae/2 Cja0(2"^x - a) for all j G Z. 

Let a; = 2-^{xo + Xi2 + 0:222 + ■■■). If j < -A^, then 2~^x G Z2, which 
implies that \2~^x — a\2 > 1 for all a G I2, a ^ 0. Thus, (f)(2^^x — a) = for all 
a G /2, a 7^ 0, and (f)(2^^x) = 1, whenever j < —N. Since, we have f{x) = Cjo 
for all j < —N. Similarly, for another x' = 2~^' (xg + a;']^2 + X222 + ■■■), we have 
f{x') = CjiQ for all j < —N'. This yields that f{x) = f{x'). Consequently, 
/(x) = C, where C is a constant. However, if C 7^ 0, then / ^ C^{Q.2)- Thus, 
C = what was to be proved. □ 

According to the above scheme, we introduce the space Wq as the orthogonal 
complement of Vq in Vi. 
Set 

(3.12) ^<«)w.^(|)_^(£_i). 

Proposition 3.3. The shift system — a),a G I2}, is an orthonormal 

basis of the space Wq. 
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Proof. Let us prove that Wo -LVq. It follows from ([31]), (KT2^ that 
(V^(°)(--a),0(--6)) = j i)^^\x-a)(t){x-h)dx 

for all a, 6 G l2- Let a ^ h. Since it is impossible a = 6+ lor6 = a + l, 
taking into account that the functions 2^/^0(2"^ ■ — c), c E I2 are mutually 
orthogonal, we conclude that [tp^^\x — a),(l){x — b)) =0. If a = 6, again due to 
the orthonormality of the system {2^/^0(2~^ ■ — c), c G I2}, taking into account 
that |, I + I G /2, we have 




Thus, + a) ± 0(- + 6) for all a, 6 G /2. 

Similarly, computing the integrals 




we establish that the system {ip^^\- — a),a E I2} is orthonormal. 
It follows from (ES]) and fl3T2D that 




If a G /2, then either a = ^b, b E I2, or a = j + -^b, b E h- Hence, 

ct>[l-a)=\{ct^{x-b)+ij^'\x-b)), bEl2, 
whenever a = ^b, and 

whenever a = | + Since {2^/^0(2"-'^ • —a) : a G I2} is a basis for Vi, we 
obtain that the system {0(- — 6), ■?/)'^°)(- — b),bE I2} is also a basis for Vi, i.e., 
the functions ip^^\- — b), b E I2, form a basis for the space Wq = Vi QVq. □ 

Thus according to Propositions l3.lt 13. 2^ 13. 3^ the collection {Vj : j G Z} is a 
MRA in £2(Q2) and the function defined by (KWf is a wavelet function. 
This MRA is a 2-adic analog of the real Haar MRA and the wavelet basis 
generated by ip^'^^ is an analog of the real Haar basis. But in contrast to the 
real setting, the refinable function (j) generating our Haar MRA is periodic with 
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the period 1 (see (13.101) ). which never holds for real refinable functions. It will 
be shown bellow that due this specific property of 0, there exist infinity many 
different orthonormal wavelet bases in the same Haar MRA (see Sec. S]). 
Due to (O, (D, the function can be rewritten in the form 

(3.13) ^(o)(x) = X2(2-M^^(|x|2). 
Thus the Haar wavelet basis is 

^(o)(x) = 2-^'/2^(°)(2^x-a) 

(3.14) =2-^'^X2{2-\Vx-a))Vt{\Vx-a\2), x G Q2, J G Z, a e h- 
Since a locally-constant function ip^^]^ [x] satisfies the relation 

(3.15) I ^fj{x)dx = 0, 

according to (12.131) . ipj^J{x) belongs to the Lyzorkin space $(Q2). 

Remark 3.1. The Haar wavelet basis f l3.14p coincides with Kozyrev's wavelet 
basis (11.21) for the case p = 2. In present paper we restrict ourself by con- 
structing the Haar wavelets only for p = 2. Since Haar refinement equation 
(13. 7p was presented for all p, a similar construction may be easily realized in 
the general case. Moreover, it is not difficult to see that Kozytev's wavelet 
function 6j{x) from (II. 2p can be expressed in terms of the refinable function 
(j){x) as 

P-i ^ ^ 

(3.16) 9k{x) = Xp{p~'kx)n{\x\p) =p-'/^J2^r<p[-x--), xeQp, 

where hr = p^^'^e^'^^^^^^ , r = 0, 1, . . . ,p — 1, = 1, 2, . . . ,p — 1. 

Remark 3.2. Because of periodicity (13.101) of the refinable function 0, one 
can use the shifts + a), a G I2, instead of ip^^\- — a), a E 

Now we show that there is another function whose shifts form an or- 
thonormal basis for Wq (different from the basis generated by tp^'^'>). Set 

,3.n, w ^ ;^(Hi) - Hi 4) - Hi - 1) - Hi 4 4)) 

and prove that the functions ^p^^\- — a) , a E I2, are mutually orthonormal. If 
a E I2, a ^ 0, |, then each of the numbers 0, ^, |, ^ -|- | differs modulo 1 from 
each of the numbers |,^ + |,| + |,^-|-| + |. Due to orthonormahty of the 
system {2^/^0(2^^x — a),a E I2} and (I3.10p . it follows that ip^^'^ is orthogonal 
to 'ip^^\- — a) whenever a E I2, a 0, |. Again due to orthonormahty of the 
system {2^/^0(2~-'^x — a), a G I2} and (I3.10p . we have 

_2-i)) = [ ij^^\x)ij^^\x - 2-^) dx 
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Thus we proved that the system {ip^^\- + a),a G I2} is orthonormaL It is not 
difficult to see that 



(3.18) 




This yields that 

^(°)(x) = i=(^a)(a;)-^«(x-i)). 

Since the system {ip^^\- — a),a G I2} is a basis for Wq, it follows that the 
system {il>''^\- — a),a G I2}, is another orthonormal basis for Wq. 

So, we showed that a wavelet basis generated by the Haar MRA is not 
unique. 

4. Description of one-dimensional 2-adic Haar bases 

4.1. Wavelet functions. Now we are going to show that there exist infinitely 
many different wavelet functions '(p^^\ s G N, in Wq generating different bases 
for C\Q2). 

In what foloows, we shall write the 2-adic number a = {clq + ai2 + ■ ■ ■ + 
as_i2^~^) G /2, flj = 0, 1, j = 0,l,...,s — 1, in the form a = where 
m = ao + ai2 + ■ ■ ■ + as-i2,'^~^. 

Since the refinable function cj) of the Haar MRA is 1-periodic (see (13.101) ). 
evidently, the wavelet function ip^ has the following property: 

(4.1) V^(°)(x± 1) = -V^(°)(a;)- 

Before we prove a general result, let us consider a simple special case. Set 

(4.2) V^«(x) = «oV^(°)(x) + «iV^(°)(x-^), «o,«ieC, 

and find all the complex numbers ao, ai for which {^^^\x — a), a G I2} is an 
orthonormal basis for Wq- 

Taking into account orthonormality of the system — a),a G I2} and 

( 14. ip . we can easily see that t/''-"^'' is orthogonal to — a) whenever a G I2, 
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a 7^ 0, |. Thus the system {ip^^\x — a), a G I2} is orthonormal if and only if 
the system consisting of the functions (14. 2p and 

(4.3) V'^'^ ~ ^) " (x) + «oV^(°) - ^) 

is orthonormal, which is equivalent to the unitary property of the matrix 

It is clear that Z) is a unitary matrix whenever |aoP + = 1- In this case 
the system — a),a E I2} is a basis for Wq because {-ip^^^- — a),a E I2} 

is a basis for Wq and we have 



So, V'*-^^ is a Haar wavelet function if and only if |aoP + = 1. In particular, 
we obtain fl3.18p for ao = «i = 

Theorem 4.1. Let s = 1,2, . . . . The function 

2'>-l , 

(4.4) ^W(3;) = ^«,^(0)(^^__), 

fc=0 

is a wavelet function for the Haar MRA if and only if 

2"-! 

(4.5) afc = 2"^(-l)'=^7,e-*-^^ A; = 0, . . . , 2^-1, 7, G C, |7,| = 1. 

r=0 

Proof Let V'^'^ be defined by (jOD. Since {V^(°)(- — a), a G is an orthonor- 
mal system (see Subsec. 13. Sp . taking into account (14.10 . we see that ip^'^^ is 
orthogonal to — a) whenever a G /2, a 7^ ^, A; = 0, 1, . . . 2** — 1. Thus the 
system {ip^'^\x — a) , a E I2} is orthonormal if and only if the system consisting 
of the functions 



(4.6) +aoV^(°) (^-^) + - ■ (^-^) 



r — 1 



, r = 0,...,2^-l. 



is orthonormal. Set 



s 



2 
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/ 
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(4.7) 



D 



\ 







^2 


Q;2«-2 


a2=-i 


— «2=-l 


ao 




a2s-3 


a2s-2 


— a2s-2 


— tt2s-l 


ao . . 


a2«-4 


028-3 




-"3 






ai 


— «! 


-«2 


-as • • 


• — a2''-i 





J 



Due to orthonormahty of {ip^^\- — a),a E I2}, the coordinates of 3^^^ form an 
orthonormal system if and only if the matrix D is unitary. 
Let u = {ao, ai, . . . , 02^-1)"^ be a vector and 





1^ 


. 


. 





-1 


\ 




1 


. 


. 










A = 





1 . 


. 



























. 


. 1 












V 


. 


. 


1 





/ 



be a 2^* X 2^* matrix. It is not difficuh to see that 
A^u = ( — a2s-r, — a2»-r+i, 



, -q;2=-i, ao, tti, . . . , a2»-r-i) ' 

T 



where r = 1, 2, . . . , 2^ — 1. Thus D = (m, Au, . . . , A^" ^u) . Hence, to describe 
all unitary matrixes D, we should find all vectors u = (ao, ai, . . . , a2»-i)"^ 
such that the system of vectors {A^u, r = 0, . . . , 2* — 1} is orthonormal. We 
have already one such a vector -uo = (1; 0, . . . , 0, 0)^ because the matrix Dq = 

(mo, Auq, . . . , A'^"~^uq)'^ is the identity matrix. Let us prove that the system 
{A^u,r = 0, . . . , 2'^ — 1} is orthonormal if and only if m = Buq, where B is 



a unitary matrix such that AB ■ 
matrix, AB = BA. Then A^u 



BA. Indeed, let u = Buq, B is a unitary 

,2" - 1. Since the 



BA'-uo, 



0,1, 



system {A^uq, r = 0, 1, 



0,1, 



1} is orthonormal and the matrix B is unitary, 
, 2^* — 1 are also orthonormal. Conversely, if 



0,1, 



2'^ — 1 is orthonormal, taking into account that 



there exists a unitary matrix B such that A^u = B{A^uq) 
Since A^' 



u 



the vectors A^u, 
the system A^u, 

{A^Uq, r = 0,l,...,2* — l}is also an orthonormal system, we conclude that 

= 0,1,. ..,^2^-1. 

A^'uq = —uq, we obtain additionally A^^u = BA^'uo- It 
follows from the above relations that {AB—BA){A'^uo) =0, r = 0, 1, . . . , 2'^— 1. 
Since the vectors A''uo, r = 0, 1, . . . , 2^* — 1 form a basis in the 2^-dimensional 
space, we conclude that AB = BA. 

Thus all unitary matrixes D are given hj D = {Buo,BAuo,...,BA'^'-^uo) , 
where is a unitary matrix such that AB = BA. It remains to describe all 
such matrixes B. It is not difficult to see that the eigenvalues of A and the 
corresponding normalized eigenvectors are respectively 
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and Vr- = {{vr)o, • • • , (fr)2»)'^5 wheie 



(4.9) 



r 

where 



0,1, 



ivr)i = 2''/\-iye~'''^', / = 0, 1, 2, . . . , 2^ - 1, 
2^ — 1. Hence the matrix A can be represented as A = CAC'^, 



A 



/ Ao . . . \ 
Ai ... 



\ ... A2=_i 
is a diagonal matrix, C = (^vq, . . . ,V2'>-i) is a unitary matrix. It follows that 
the matrix B = CBC~^ is unitary if and only if B is unitary. On the other 
hand, AB = BA if and only if AB = BA. Moreover, since according to (14.81) . 
Afc 7^ A/ whenever k I, all unitary matrix B such that AB = BA, are given 
by 

/ 70 ... \ 
7i ... 

B = 

V ... 72.-1 / 
where 7^ G C, = 1. Hence all unitary matrix B such that AB = BA, are 
given hy B = CBC^^. Using ( 14. 9p . one can calculate 

2"-! 

ttfc = (-BMo)fc = {CBC~^UQ)k = ^ 'JriVr)kiVr)o 



r=0 



2"-! 



k = 0,l. 



r=0 



where 7^ G C, |7fc| = 1. 

It remains to prove that the {ip^^\- — a), a G I2} is a basis for VFq whenever 
ijj^^^ is defined by ( Hill) . fH3|) . Since {^/;(°)(- - a), a G 12} is a basis for Wq, 
it suffices to check that any function ■?/'(°)(- — c), c G I2, can be decomposed 
with respect to the functions ip^^\- — a), a E l2- Any c G /2, c 7^ 0, can be 
represented in the form c = ^ + b, where r = 0, 1, . . . , 2'^ — 1, \b\2 > 2*"'"^. 
Taking into account that = i.e.. 



2^-1 



^■•°'(--^) = E'3"*.'"(x-| 



r = 0,l,...,2" 



we have 



(0) 



X — c 



fc=0 



^(0) ( a; _ _ _ 5 



2=-l 



fc=0 



2^ 



and ;| + 6g/2, A; = 0,l,...,2^ 



1. 



□ 
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Finally note that, due to fl27[3|) . /^^ ^^'^(2^ 

x—a) dx = 0, j G Z, a G I2, which 
yields that any function ■ip^'^\2^ ■ —a) belongs to the Lizorkin space ^(Qg). 

4.2. Real wavelet functions. Using formulas ( 14. Sp . one can extract all real 
wavelet functions (14.41) . 

Let s = 1. According to (Oil . 

(4.10) V^W(x) = cos^^(°)(x) +sin0^(°)(^x- ^) 

is the real wavelet function. 

Let s = 2. Set 7^ = e*^^ r = 0, 1, . . . , 2" - 1. It follows from (gSD that the 
wavelet function ijj^^^ is real if and only if 

sin 61 + sin 62 + sin 63 + sin 64 = 0, 

cos 9i — cos 02 + cos 6^3 — cos 04 = 0, 

sin 61 — sin 62 — sin 63 + sin 64 = 
cos 6*1 + cos 6*2 — cos 6*3 — cos 6*4, 

sin 61 — sin 62 — sin ^3 + sin 64 = 
— (cos 6*1 + cos 6^2 — cos6'3 — cos6'4). 

These relations are equivalent to the system 

sin^i = — sin^^4, cos^^i = cos 6*4, 
sin 6*2 = —sin 6*3, cos 6^2 = cos 6*3. 

Thus the real wavelet functions (14. 4p with s = 2 are given by 
V;«(x) = ^(cos^i + cos^2)V^^°^(a;) 

+ ^^ (cos 61 — COS 6*2 + sin^^i + sin 6*2)'?/'^°^ {x — — j 

+^(sin^i - sine2)^^°^ [x - 

(4.11) + (cos 6^1 - cos 6^2 - sin 6^1 - sin6'2)?/'^°^ fx - ^ - - 

2v2 V 2^ 2. 

In particular, for the special cases 9i = 62 = 0, 61 = —62 = 6, 61 = 62 + ^ = 6, 
we obtain respectively the following one-parameter families of real wavelet 
functions 

^(1) (x) = cos (x) + sin ^V^W (a; - 1) , 
(x) = cos (x) + ^ sin ^V^(o) (^x - ^) 

(4.12) _^sine^(o)(^a;_^_ 1), 

tjj(^\x) = i(cos^-sin^)^/'(°)(x) +^(cos^ + sin^)^/'(o)(^x-^) 

-i(cos^-sin^)?/^(o) C 
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5. Description of multidimensional 2-adic Haar bases 

5.1. p-Adic separable multidimensional MRA. Here we describe multi- 
dimensional wavelet bases constructed by means of a tensor product of one- 
dimensional MRAs. This standard approach for construction of multivariate 
wavelets was suggested by Y. Meyer [29] (see, e.g., [30i §2.1]). 

Let {V}%^^ , z/ = 1, . . . , n, be one-dimensional MRAs (see Subsec. 13.11) . 
We introduce subspaces Vj, j E Z, of £^(Qp) by 

n 

(5.1) V, = (g) Vl"^ = span{F = /i ® ■ ■ ■ ® /„, /. G vj"^}. 

Let be a refinable function of z/-th MRA {vf^}j. Set 

(5.2) $ = 0W®---®0("). 

Since the system {0'^ '')(■— a) ja^^g/^ is an orthonormal basis for V^^^ (axiom (e) 
of Definition 13.11) for any = 1, . . . , it is clear that 

Vo = span{<l>(- - a) : a = (ai, . . . , a„) e /^}, 

where 1^ = Ip x ■■■ x Ip is the direct product of n sets Ip, and the system 
$(■ — a), a G Ip, is an orthonormal basis for Vq. It follows from Definition 
fl5T|) and axiom (d) of Definition O that f e Vq if and only if /(2-J -) e Vj for 
all j G Z. Since axiom (a) from Definition 13.11 holds for any one-dimensional 

MRA {Vj^^jj, it is easy to see that $(2"^' ■ -a) G for any a E I^. Thus, 
Vj C It is not difficult to check that the axioms of completeness and 

separability for the spaces Vj hold. Thus we have the following statement. 

Theorem 5.1. Let {vj^^jj^z, u = I, . . . ,n, be KMAs m ^^(Qp). Then the 
subspaces Vj of C'^{Qp) defined by (15.11) satisfy the following properties: 

(a) Vj C Vj+i for all j G Z; 

(b) Uj^zVj IS dense m C^iQp); 

(c) n.^zV, = {0}; 

(d) /(■) G V, ^ f{p-'-) G V^,+i for all] G Z; 

(e) the system {$(x — a), a G Ip}, is an orthonormal basis for Vq, where 
$ G Vo zs defined by (1^ . 

Similarly to Definition 13.11 the collection of spaces Vj, j G Z, which satisfies 
conditions (a)-(e) of Theorem 15. II is called a multiresolution analysis in £^(Qp), 
the function $ from axiom (e) is called refinable. 

Next, following to the standard scheme (see, for example, pHi §2.1]), we 
define the wavelet spaces Wj as the orthogonal complement of Vj in V^+i, i.e., 

Wj = v,+ieVj, jGZ. 

Since 

n n 
v=l u=l 
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eC{l, .■■,"}, e^0 '^C'^ M^e 

So, the space ly^- is a direct sum of 2" — 1 subspaces Wj,e, e C {1, . . . , ra}, e 7^ 0. 
Let t/''-'^^ be a wavelet function, i.e. a function whose shifts (with respect to 
a e Ip) form an orthonormal basis for Wq''\ It is clear that the shifts (with 
respect to a G ) of the function 

(5.3) *e= ((g)^^^O((g)0^'^)), ec{l,...,n}, e ^ 0, 

form an orthonormal basis for Wo,e- So, we have 

iez jez ec{i, .■■,"}, 67^0 

and the functions j9~"-'/^\E'e(p^ ■ +a), e C {1, . . . , n}, e 7^ 0, j G Z, a G , form 
an orthonormal basis for £^(Qp). 

5.2. Construction of multidimensional 2-adic Haar MRA. Let us apply 
the above construction taking the 2-adic Haar MRA as u-th one-dimensional 
multiresolution analysis {Vj'^^}j(zz, v = 1, . . . , n. 

To construct multivariate wavelet functions (15.31) . we choose ■^/^'-^^ as a wavelet 
function for each one- dimensional MRA. Thus we have the following 2" — 1 
multidimensional wavelet functions 



= ^(°)(a;i)0(a;2) ■ ■ ■ 0(x„_i)0(x„), 

^J°\ = 0(X1)0(X2) ■ • ■ 0(X„„1)7/;(°)(X„). 

Let e C {l,...,n}, e 7^ 0. Denote by ke = {{ke)i, ■ ■ ■ , {ke)n) the vector 
whose coordinates are given by 

(fc,).= |j^ l^l^^ u=l,...,n. 

Since (f){x,y) = ^l{\xu\2) and ^/^(^^(x;,) = X2{'2~^x,y)^l{\x,y\2), Xy G Qp, z/ = 
1, 2, . . . , n, (see (jSlED and fl313D ). it follows from (Q, ([31]), that the wavelet 
function '^f^ can be rewritten as 

(5.4) ^f(x) = X2(2-^A;e-x)fi(|x|2), x = (xi, . . . , x,) G Q^. 
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According to the above consideration, we have we the following statement 
Theorem 5.2. The system of functions 

(5.5) =2-"^/2^2(2"^A;e- (2^a;-a))r](|2^a;-a|2), x G Q^, 

e C {1, . . . e 7^ 0, j E^L, a E is an orthonormal basis for i2^(Q2). 

Now we construct multidimensional wavelet bases using different one- dimen- 
sional Haar wavelet bases (see Subsec. I4?TI) . Namely, we apply the construction 
of Subsec. 15.11 taking again the Haar MRA as one-dimensional multires- 
olution analysis {vj'^^}j^z, ^ = 1, . . . , ri, and choosing wavelet functions ip'^^"'^ 
for construction of multivariate wavelet functions ( 15. 3p . 

Let s = (si, . . . , s„), where G No, v = 1, 2, . . . , n. We have the following 
2" — 1 wavelet functions 

*{l!...,n}(^) = i^^''\^l)i^^''\^2) ■ ■ ■^(^-^H^n-OV'^'-na^n), 
^g,,,^„_^^(x) = ^(^^)(X1)^(^^)(X2) ■ ■ ■^(^"-^Ha;„-l)0(a;n)^ 



^^^Ux) = ^IJ^''\xi)<P{x2) ■ ■ ■ 0(a;„_i)0(x„), 



^J;\(X) = 0(X1)0(X2) ■ ■ ■ <P{Xn-M^'-\Xn). 

Set al = ar, where is given by (14. 5p . r = 0, 1, . . . , 2* — 1, and «[} = 1, 
a° = ■■■a2s_i = 0. Since 0(xy) = i7(|xjy|2), Xi, G Q2, and ■?/;('^'^) is given by 

(I44|) . (I43|) . u = l,2,...,n, due to (12151) . the wavelet functions ^i''^ can be 
rewritten as 

v^(^)(x) = 

(5.6) 

2=1-1 2="-l 

= E---E«S^^^---<^^"*i°^(--(^(^e)i,...,^(M.)), xGQ^, 

n=0 r„=0 

where \Efe°'' is defined by ( 15. 4p . e C {1, . . . , n}, e 7^ 0. 

According to the above consideration, we have the following statement. 

Theorem 5.3. The system of functions 
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(5.7) = E E ^i''''---^il''"'^'ei2^--a-{^iK)u..^^ 

ri=0 r„=0 

X e Q2 J e C {1, . . . , n}, e 7^ 0, j E Tj, a E I2 forms an orthonormal basis for 



Since for a locally-constant function \l'gy-^(a;) the relation J„„ \I^g'^]^(x) 



x 



holds, in view of fl2:T3|) . G $ 



6. P-AdIC WAVELETS AS EIGENFUNCTIONS OF PSEUDO-DIFFERENTIAL 

OPERATORS 

6.1. Pseudo-differential operators. S. V. Kozyrev [20| proved that one- 
dimensional p-adic wavelets (11. 2p are eigenf unctions of the Vladimirov frac- 
tional operator a > 0. In fact, this statement holds for all a G C due 
to (I3.15p . A criterion for pseudo-differential operators (12.151) to have wavelets 
(II. 2p and (11.31) as eigenf unctions was found in [5], [IB]. 

Now we consider a similar problem for 2-adic wavelets (15. 5p and (15. 7p . 

Theorem 6.1. Let A be a pseudo- differential operator (I2.15P with a symbol 
A G ^(Q^' \ {0}), e C {1, . . . , ra}, e 7^ 0, j eZ, ael^. Then the function 

^g,(x) = 2-^^/\2{2^'h ■ - a))n{\2^x - a^), x G 
is an eigenfunction of A if and only if 

(6.1) A{V{-2'^K + r])) = A{-V-^K), Vr/GZ^. 

The corresponding eigenvalue is X = A[ — 2^~^k(,), i.e., 



r2, 



A^Z = Ai-2^~X)¥:j^. 
Proof Combining ^Alh . fl2T2|) . fl23|) with ([53]), we obtain 
i^[^Sj(0 = 2-^/'F[¥^\2^x - a)](0 = 2-^/\2{2-'a ■ ^)F[¥:>\x)]{2-^0 

n 

= 2"^-/2^2(2-%-e)i^[n^2(2-^(M.a;.)fi(|x.|2)j(2-^0 



11 

2-^/\2{2-'a ■ ^)l[F\n{\x,\2)] {2-\K), + 2'^^ 



(6.2) = 2"^/2^2(2-% ■ 0^(|2"'fce + 2-^'e|2). 

It is clear that n[\2-^{ke)y + ^^.b) ^ only if 



-2 H^e)i. + Vu, where 



77^ G Z2, z/ = 1, 2, . . . , n. This yields ^ = -2^1/^^ + r^, G Z;j. 
If condition (16. ip is satisfied, then, using (I2.15p . (16. 2p . we have 



(6.3) 



2".72^"i [^(e)x2(2-^a • 0fi(|2-'fce + 2-^i\,)] 



X] 
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Making the change of variable ^ 2^{—2^^kf. + r]) and using fl2.12p . we obtain 
Avl>g„(x) = 2-"^/2y ;^2(-(2^x-a)-(-2-iA;,+r7))^(2-'(-2-iA;,+r7))r](|r7|2)c?"r/ 

= 2-'^^/'A{-2^-'j)x2{2-% ■ {2^x - a)) [ X2(-(2^x - a) ■ v) 

= A{-2^-'j)^'^l{x). 

Consequently, A¥^j^{x) = A^J^^al^^)' where A = Ai-2^~^K). 

Conversely, if ^^E^g^ja = A\l/g°j^, A G C, taking the Fourier transform from 
both the left and the right hand sides of this identity and using (I6.2p . (I6.3p . 
we have 

(6.4) (^(0-A)x2(2-%-e)^^(|2-^A;e + 2-^e|2) =0, ^ e Ql 

If now G Z^, ^ = 2-'(-2-iA;e + v), then 2-iA;e + 2'^^ = r/ G Z^. Since 
fi(|2-iA;e+2-^'^|2) ^ 0andx2(2-^a-^) ^ 0, it follows from ([63D that A = ^(0- 
Thus A = ^(2^(-2-iA;e + v)) for any rj G Z^. In particular, A = A{-2^~^j), 
and, consequently, (16.11) holds. □ 

In particular, we have the following statement. 

Corollary 6.1. Let e C {1, . . . , n}, e 7^ 0, j G Z, a G /2 ■ Then the function 
^e°]a ^-^ '^'^ eigenfunction of the fractional operator (12.231) . The corresponding 
eigenvalue is X = 2"^^~^\ i.e., 

Proof. The symbol A{C) = |^|2 of the fractional operator satisfies condition 
(!6T|): 

^(2^(-2-^A;e + V)) = |2^(-2-iA;e + = 2"'"( max | - 2-i(fce). + VA2T 
= 2-(i-^)( max |2^(-2-i(A;e),)L)" = ^(-2^-ifce) 

^ l<u<n ' 

for all rj G Z2, e C {!,...,«}, e 7^ 0. Here we take into account that 
{ke)u = 0, 1; = 1, 2, . . . , n; (fce)i + ■ ■ ■ + (/Ce)n 7^ 0. Thus, by Theorem 16.11 
\l/g°jjj is an eigenfunction and the corresponding eigenvalue is A = 2°^^^'^^ □ 

Similarly to Theorem 16.11 and Corollary 16. H using representation (15. 7p . it is 
not difficult to prove the following statements. 

Theorem 6.2. Let A be a pseudo- differential operator (I2.15P with a symbol 
A{0 e ^(Q^ \ {0}), s = {si,...,sn), where G No, e C {1, . . . , n}, e ^ 0, 
j G Z, a G ■ Then the function ^I/g^ja is an eigenfunction of A if and only if 
(16. ip holds. The corresponding eigenvalue is X = A[ — 2^~^ke), i.e., 

A^% = A{-2^-X)¥:l. 
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Corollary 6.2. Let s = (si, . . . , where s^, G Nq, e C {1, . . . , n}, e 7^ 0, 

is] 

j G Z, a G Then the function '^),.ja is an eigenfunction of the fractional 
operator fl2.23p . The corresponding eigenvalue is \ = 2°'^^^^\ i.e., 
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